Dipole Oscillations of a Bose-Einstein Condensate in Presence of Defects and Disorder 
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We consider dipole oscillations of a trapped dilute Bose-Einstein condensate in the presence of 
a scattering potential consisting either in a localized defect or in an extended disordered potential. 
In both cases the breaking of superfluidity and the damping of the oscillations are shown to be 
related to the appearance of a nonlinear dissipative flow. At supersonic velocities the flow becomes 
asymptotically dissipationless. 
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One of the most spectacular consequences of phase co- 
herence and interactions in condensed matter is super- 
fluidity, a direct manifestation of which is the capacity 
of a fluid to move without dissipation. According to 
the standard Landau criterion, the superfluidity (SF) of 
a uniform flow of, e.g., liquid He 4 , or a Bose Einstein 
condensate (BEC) will be broken if an obstacle moves 
through the fluid with a speed higher than a critical 
velocity v\ = min {E(p)/p} where E(p) is the disper- 
sion relation of elementary excitations of momentum p. 
Though this property has been explicitly checked in He 4 
[H and in a BEC flow Q in presence of small impurities, 
experiments in superfluid 4 He (see, e.g., [3]) and more 
recently in BEC [3, H[ have shown that the critical ve- 
locity for breaking SF is generically lower than vf?, due 
to phase slips induced by vortex (or soliton) emission, as 
originally proposed by Feynman [6[ . 

Collective oscillations of BEC confined by harmonic 
traps offer new opportunities to explore the central ques- 
tion of the breaking of SF and of the origin of drag and 
dissipation in quantum liquids and gases. In a recent se- 
ries of experiments, damping of the oscillations (such as 
dipole or quadrupole) in the presence of a single local- 
ized scatterer [7[, and disordered [8, 9] or quasiperiodic 
[Tot superimposed potentials has been used to investigate 
different dynamical regimes, including the possibility of 
a Bose glass, (Anderson) localization or other possible 
phases. These investigations have clearly shown the ex- 
perimental relevance of analyzing transport properties of 
BEC via the damping of collective excitations. However 
the connexion of the damping with localization proper- 
ties still remains to be clarified. 

Our purpose here is to provide a global analysis of the 
phase diagram related to the damping of dipole oscilla- 
tions in the presence of a single localized scatterer or a 
disordered potential. We consider the regime where the 
experiments have been realized up to now, i.e., a quasi- 
1D geometry where the chemical potential is larger than 
the typical amplitude of the perturbing potential. The 
reason why we treat together the localized defect and 
the random potential is that, qualitatively, several of the 
main features of the dynamics are contained in the former 




FIG. 1: (color online) Dynamical regimes for dipole oscilla- 
tions in presence of a localized Gaussian defect. The plot 
represents the fluidity factor 7 (see the text) computed after 
a time tf = 25 x 2tv/oo x . The yellow (light gray) region corre- 
spond to zero damping (7 = 1). The dashed lines are analytic 
determinations of the frontiers between the different regimes 
(see text). 



case. Its analysis therefore facilitates the comprehension 
of the latter, and stresses the generic aspects, leading to 
a unified picture of dissipation. We find that in both 
cases there exist SF undamped oscillations at small am- 
plitudes. As the amplitude of oscillation (or the typical 
size of the perturbation) increases, the system enters a 
dissipative regime where solitons and phonon-like excita- 
tions are emitted (a regime recently studied experimen- 
tally in Ref.[5[ for a moving localized obstacle). Though 
the system preserves its phase coherence, the collectivity 
of the center of mass motion, and therefore the amplitude 
of the dipole oscillations, diminishes [8|, [9|. In the case of 
a superimposed disordered potential, this dissipative or 
resistive phase, where non-linearities of the system play a 
crucial role, has no relation with (Anderson) localization. 

The system considered is a weakly interacting BEC 
confined in a cylindrically symmetric 3D harmonic po- 
tential m{iJ\r\ + c^£ 2 )/2 in presence of an additional 
potential U(x). In the limit of a highly anisotropic trap, 
cu± ^> lj x , the transverse confinement is such that the 
quasi- ID regime can be reached. It is important to note 
that for moderate U (x) (even a disordered one) the phase 
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coherence of the system is preserved as demonstrated in 
Refs. [|, [Tl|. The system is thus accurately described 
by a ID order parameter ip(x,t), depending on a single 
spatial coordinate x along the axial direction of the trap. 
ip(x,i) obeys the nonlinear Schrodinger equation [12L [l3(] 
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Here, n{x,t) = \i(j(x,t)\ 2 is the condensate density per 
unit of longitudinal length and a > is the 3D s- 
wave scattering length. In the low density regime (LDR, 
an < 1) the density profile in the transverse direction is 
Gaussian-shaped and v = 1, whereas v = 1/2 in the op- 
posite high density regime (HDR, an ^> 1) where the 
Thomas-Fermi approximation holds for the transverse 
degree of freedom. Note that Eq. (pQ) does not account 
for transverse excitations which may be relevant in the 
HDR. We checked that it nonetheless gives an excellent 
account of the experimental result on dipole oscillations 
of the Florence and Rice groups [8|, [9| performed in the 
HDR. 

After preparing the condensate in the ground state of 
the harmonic trap (with density tiq(x), chemical poten- 
tial a and, at the center of the cloud, speed of sound c), 
dipole oscillations were excited by a sudden displacement 
do of the harmonic potential. For U(x) =0 the center of 
mass oscillates freely with frequency w x , and acquires a 
velocity v = co x do when passing through the origin. The 
time evolution of the density reads n(x,t) = no(x — X t ) 
where X t = docos(u x t) is the position of the center 
of mass. For a finite U(x), which is turned on simul- 
taneously with the sudden displacement of the trap, 
X t = f R dxxn(x,t) is computed numerically up to 
a time tf chosen such that X t>tf assumes an oscillatory 
pattern of roughly constant amplitude which we denote 
df. In order to measure the damping of the dipole oscil- 
lations we define a fluidity factor 7 = df/do (7 = 1 in 
the absence of damping and 7 — > for strong damping). 

Localized defect. We start by considering a Gaussian- 
shaped defect U(x) — Uo exp(— x 2 /2<j 2 ). The fluidity 
factor 7 is plotted in Fig. [J as a function of the nor- 
malized defect strength Uo/u and velocity v/c = uo x do/c. 
The numerical calculations were performed for a BEC 
in the LDR with chemical potential \i — 40 hw x and as- 
pect ratio oo x /co± selected to have a Thomas- Fermi like 
density-profile along the axial direction. In this case 
no(x) = M v G(L - \x\) [L 2 - x 2 \ x ' v where the factor M v 
normalizes the density to the number of atoms N, Q(x) 
is the Heavyside step function and L = ^/2ji/mou 2 is half 
the longitudinal size of the condensate. The parameters 
are TV = 1.5 x 10 4 87 Rb atoms, uj x = 2tt x 9 s" 1 = uj±/10 
and a = £ = 0.28 /im, where £ = h/^2mu is the healing 
length at the center of the condensate. The qualitative 
structure of Fig. [1] is generic and does not depend on the 
specific values of a and /i, and is also observed in the 
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FIG. 2: (color online) Density profile after a time t — 3/4 x 
2ty/co x for Uo//n = 0.24 at different initial velocities: a) v/c = 
0.1, b) v/c = 0.67, c) v/c = 1.2, d) v/c = 2.5. The confining 
potential is represented as a full (red) curve. Insets blow up 
the density around the defect. (Other parameters as in Fig.[T]) 



HDR [y = 1/2). 

In the deep subsonic limit v/c <C 1, the Gaussian scat- 
terer induces no observable damping of the dipole oscil- 
lations. Numerical results show that the oscillating con- 
densate is only locally perturbed in the vicinity of the 
defect: a dip or a peak appear in the condensate den- 
sity for L^o > and Uo < 0, respectively (see Fig. [2fa)). 
These are characteristic features of a superfluid flow, with 
no energy dissipation, nor drag exerted [14], and with no 
damping of the oscillations (perfect transmission through 
the scatterer potential). In this regime, a perturbative 
treatment of equation (pQ) and a local density approxima- 
tion yield a condensate density of the form 



n(x, t) = no(x — X t ) [1 + 5n(x, t)] 
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where k = ^|cq — J*Q 2 | 1//2 , Co being the unperturbed 
local sound velocity: mc\ — 2vhw±_[ano{x — X t )) v . The 
accuracy of this approximation is shown in Fig. EJa); it 
is well justified if mUocr/(h 2 hz) <C 1 when kg <C 1, or 
mUo/{hn) 2 < 1 when kg > 1 [13], and if g <C L. 

For weak defect potentials, \Uo\ <C /i, and if the TF 
size L is large compared to the dipole oscillation ampli- 
tude, the center of mass position X t can be computed 
analytically. To lowest order, the solution of the small- 



amplitude linearization yields X t = do cos[(uj x 
where the defect-induced frequency shift reads 
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defect, in excellent agreement with our numerical results. 
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Hence the analytical evaluations of the density profile ([2]) 
and of the center of mass motion confirm the superfluid 
behavior of the oscillations in the deep subsonic regime. 

The situation changes as the velocity increases at fixed 
Uo/fi or as Uo/fi is increased at constant velocity. In the 
former case, at some critical velocity v c (< c), that de- 
pends on the strength of the defect potential, the system 
looses SF, damping is observed and the fluidity factor di- 
minishes. In the spirit of Landau's criterion we identify 
the border between the SF and this "dissipative region" 
as the locus of points where the maximum local con- 
densate velocity v(x,i) equals the local speed of sound 
c(x,t) = (2uhujj_/m) 1 ^ 2 [an(x,t)] u ^ 2 [TBI] . The former can 
be computed from mass conservation and Eq. ([2|). For 
an impurity localized at x = this yields 

^ = [1 + Suc) 1 ' 2 ^'" if U > , - = 1 if U < 0. (5) 

c c 

where Sn c is the factor 5n(x, i) [Eq. (j3])] evaluated when 
x = X t = 0. Fig. [1] shows that these estimates coincide 
well with the numerical findings. 

When the system enters the dissipative regime, one or 
a few gray solitons detach from the defect during the first 
oscillations, as well as some pho non-like excitations (see 
Fig. [2Kb)). As time goes on, the interactions of the soli- 
tons among them, with the defect and with the phonon- 
like excitations produce time-dependent fluctuations of 
the shape. During this process the condensate does not 
loose phase coherence, but the center of mass motion 
looses collectivity, part of the kinetic collective energy 
being transformed into density fluctuations. The damp- 
ing process continues until the center of mass velocity 
becomes comparable with the critical velocity. There- 
upon, though presenting local density fluctuations, the 
amplitude of the oscillations remain constant in time. 
Deeper in the dissipative regime, an increased emission 
of gray solitons and phonon-like excitations is observed, 
leading to a massive distortion of the initial condensate 
profile (see Fig. Etc)). The BEC enters a strongly irreg- 
ular time-dependent regime, the collectivity of the dipole 
motion is totally lost, and the damping increases drasti- 
cally. 

Finally, at sufficiently high supersonic velocities, a dif- 
ferent phase is reached where the damping tends again to 
zero (7 — > 1). In this regime, that we denote as "quasi- 
ideal" , the kinetic energy of the condensate is large com- 
pared to the strength of the external potential, and inter- 
actions tend to be negligible. We find, in agreement with 
previous theoretical studies [HI, [3| ? a strong suppression 
of dissipation as the velocity increases. The condensate 
density is again only locally distorted in the vicinity of 
the defect (cf Fig. Etd)). This local distortion is very 
well described by applying the combination of perturba- 
tive and local density approach already used in the SF 
regime. The density is of the same form as in Eq. ([2]) 
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FIG. 3: (color online) Fluidity factor 7 of dipole oscillations 
in presence of a speckle potential. 
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where U is the Fourier transform of U. This form indeed 
corresponds to almost perfect transmission, with small 
reflexion on the defect (see the inset of Fig. Etd)), the 
amount of which decreases at large velocity (k — > 00). 

The frontier between the dissipative and the quasi- 
ideal region can be estimated by studying the related 
problem of a homogeneous condensate flowing through a 
barrier potential [16]. In this simplified configuration it 
is possible to determine analytically the velocity at which 
the system undergoes a transition from a local perturba- 
tion to an irregular fluctuating density profile [l3j . We 
find that this estimate fits very well the numerically de- 
termined supersonic frontier (see Fig. [1]). This stresses 
the qualitative similarities between Fig. [1] and the phase 
diagram obtained for a defect moving through a homo- 
geneous fluid [HI, Q3] • Interestingly, the existence of the 
three regions (SF, nonlinear dissipative, and quasi-ideal 
weakly-damped) was quite recently observed experimen- 
tally for a localized defect in Ref. [5[. 

Disordered potential. Keeping the same parameters as 
in Fig. [U we now replace the single localized impurity by 
a disordered potential and compute, as before, the fluid- 
ity factor as a function of velocity and of the intensity of 
the (now random) potential. Fig. [3] corresponds to the 
case where U(x) is an optical speckle potential of mean 
value U with a correlation length l c such that L/l c = 30, 
typical of experimental configurations [17[. The picture 
is generic and does not depend of the details of the dis- 
order. The main result that emerges from the compar- 
ison of Fig. [D and Fig. [3] is that, in the weak-disorder 
limit that we explore (U / ji <C 1), the global properties of 
the damping phase diagram are qualitatively similar in 
both cases. The same three phases observed for the lo- 
calized defect are again present. However, their relative 
importance is quite different. One observes a consider- 
able shrinking of the SF and quasi-ideal weakly-damped 
regions, compared to the nonlinear dissipative one. 
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FIG. 4: (color online) Relative standard deviation 
(Sou 2 ) 1 / 2 as function of average speckle intensity. The 
analytical results correspond to Eq. J7}. The insets compare 
numerical and analytical frequency shifts for different realiza- 
tion of disorder for U / ji — 0.035 in the HDR and LDR. 



In the presence of disorder, important experimental ef- 
forts have been undertaken for understanding the possi- 
ble connection of damping of dipole oscillations with An- 
derson localization of matter waves. Our analysis shows 
that damping occurs in the dissipative phase. However, 
the dissipative mechanism observed in this phase is -as 
in the case of a localized defect- connected to the loss 
of collectivity due to the emission of solitons and lin- 
ear excitations, and not to a localization phenomenon. 
From what is known in the case of an a homogeneous 
flow, genuine Anderson localization might only occur in 
the deep supersonic regime v/c ^> 1 [18]. Without go- 
ing into a detailed analysis of the experimental results 
[1, [H , let us mention that the displacement of the har- 
monic potential (around 700 jam) used at Rice University 
corresponds to v/c ~ 2.8, while the lowest speckle height 
considered is U/fi « 0.04, which locates the system in the 
dissipative phase (under the experimental conditions, we 
find that at v/c ~ 2.8 the dissipative phase border is 
at U/a « 0.008). Similarly, from the data published in 
Ref. [8] our analysis shows that the experiments at Flo- 
rence were also performed in the supersonic dissipative 
region. These simple remarks explain the experimentally 
observed damping, and locate the experimental configu- 
rations in a regime where the effects of Anderson local- 
ization are at best indirect. 

In the subsonic SF regime it is possible to compute the 
frequency shift due to the disordered potential similarly 
as in the one-peak case. The approach is found to be 
equivalent to the sum rule approach developed in [1, Qjl . 
One can also derive a simple relation for the variance of 
the frequency shift: 
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where (...) denotes ensemble average. For small U//jl 
(when the lowest order approximation holds) Eq. (0) is 



found to be in very good agreement with numerical inte- 
gration of the GP equation, as seen in Fig. [4| 

To conclude, we have presented a comprehensive pic- 
ture of the damping properties of dipole oscillations of 
BEC in the presence of a scattering potential. Strong 
analogies are stressed between different types of poten- 
tials. Three different phases are shown to exist: super- 
fluid (v/c < 1), nonlinear dissipative (v/c ~ 1) and quasi- 
ideal (v/c > 1). The mechanism that breaks SF and 
leads to damped oscillations is shown to correspond to 
a generic onset of dissipation which, in the presence of 
disorder, is unrelated to localization properties. As the 
strength of the disorder potential increases the nonlinear 
dissipative phase occupies most of the phase diagram. 
Our findings allow to give a simple interpretation of ex- 
perimental results. 
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